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I. INTRODUCTION 

Since LEPS collaboration reported the 0 + pentaquark 
Q, many subsequent experiments claimed the confirma¬ 
tion of this exotic baryon with S=+l while many other 
groups didn’t observe it in their search. Experiments 
with both positive and negative results are reviewed in 
Refs. J2], 3j- At present, neither the existence nor the 
non-existence of 0 + is established, which can only be 
settled by the new high-statistical experiments. 

If 0 + really exists, its low-lying mass, extremely nar¬ 
row width and weird production mechanism pose a se¬ 
rious challenge to theorists. There have been over four 
hundred theoretical papers addressing these issues. The 
early theoretical development can be found in a recent 
review Q. The angular moment of 0 + is assumed to be 
J = i in order to render 0 + low-lying in most of these 
efforts. Moreover, some models prefer the positive parity 
to ensure a narrow 0 + . 

However, the possibility of J = | is not excluded. For 
example, J p = ^ + pentaquarks are always accompanied 
by J p = | + partners Q in Jaffe and Wilcezk’s diquark 
model [lj. Their magnetic moments were calculated in 
jj] while their radiative decays and photoproduction were 
studied in [|j. 

If 0 + carries J p = | , it decays into nucleons and 
kaons via D-wave. The phase space suppression factor 

is ~ ~ 10 -5 where px is the decay momentum 

of the kaon. The decay width could be well below ten 
MeV even if the coupling constant Gqkn is big due to 
the absence of the orbital excitation inside 0 + and 0’s 
strong overlap with KN. 

There have been a few theoretical papers on the possi¬ 
ble J = | pentaquarks using different models. Page and 

Robert suggested I = 2, J p = \ , | for 0 + to resolve 
the narrow width puzzle fg]. Jaffe and Wilczek discussed 
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the J p = | assignment for S pentaquark 0- The 

mass spectrum of J p = § pentaquarks were studied 
with the perturbative chiral quark model a Takeuchi 
and Shimizu suggested the observed 0 resonance as a 
I = 0, J p = | NK* bound state using the quark model 
pH ]. With the flux tube model, Kanada-Enyo et al. stud¬ 
ied the mass and decay width of the I = 1, J p = | 
pentaquark jp . Huang et al. proposed 0 + as a molec¬ 
ular state of NKn with 1 = 1, J p = § using the chiral 
SU(3) quark model 0. The phenomenology of J = | 
pentaquarks such as the mixing scheme and mass pattern 
was discussed in 0. Pentaquark states with J=3/2 and 
1=0,1 were studied using currents composed of one scalar 
diquark and one vector diquark .16) ■ 

We shall employ QCD sum rules (QSR) to explore the 
possible existence of the J = | ,S = +1 pentaquark 
states with the isospin I = 0,1,2. QSR formalism was 
first employed to study the pentaquark mass with differ¬ 
ent isospin in Ref. [13 ■ Up to now there have been more 
than ten papers on pentaquarks within this framework 

ElillialiSElElEllilElllliliiil. 

In practice there are two forms of QCD sum rules. 
The first one is the conventional Laplace sum rules intro¬ 
duced originally by the inventors of this formalism ^p]. 
The other one is the finite energy sum rule (FESR) |3ll |. 
Their difference lies in the weight function. The right- 
hand-side (RHS) of the traditional sum rules deals with 
fs° p{s)e~ s / MB ds, where p(s) is the spectral density 
including the nonperturbative power corrections arising 
from various condensates. Mb is the Borel parameter. 
Smin is the starting point of the integral, which is zero 
for massless quarks. In the analysis of the sum rules, the 
quark-hadron duality is always invoked. Starting from 
So, the physical spectral density, which arises from the 
higher resonances and continuum is always replaced by 
the perturbative one. Hence so is called the threshold 
parameter and is typically around the radial excitation 
mass. 

For the FESR approach, the exponential weight func¬ 
tion e~ s ! M B is replaced by s n in the numerical analysis. 
For the conventional ground-state hadrons such as the 
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rho and nucleon, both the Laplace QSR and FESR yield 
almost the same numerical results for the hadron mass, 
thanks to (1) the good convergence of the operator prod¬ 
uct expansion; and (2) the useful experimental guidance 
on the threshold parameter so- The reason is simple: the 
rough value of the radial excitation is more or less known 
exp eriment ally. 

In the present case, the situation is more tricky. Even 
the existence of the lowest 0 + pentaquark has not been 
established, let alone its radial excitation. On the other 
hand, the spectral density p(s) ~ s m with m > 5, 
which causes strong dependence on the continuum or the 
threshold parameter s 0 - Hence, compared to the Laplace 
sum rule with the double parameters (so, Mg), the single¬ 
parameter FESR may have some advantage as recently 
noted in [gjj. With FESR, one can study the dependence 
of me on so in the working region. 

On the other hand, the disadvantage of FESR is that 
its weight function s n enhances the continuum part even 
more than the weight function e~ s ! M s j n the Laplace 
sum rule. Hence some uncertainty is connected to the 
continuum threshold sq- Especially one must make sure 
that only the lowest pole contributes to the FESR be¬ 
low so- Otherwise the result will be misleading, which 
will be shown explicitly in our numerical analysis. To be 
more specific, a naive stability region in so is no guaran¬ 
tee of a physically reasonable value for sq- For example, 
the FESR with an extracted threshold So ~ 20 GeV 2 is 
certainly irrelevant for the 0 + pentaquark around 1.53 
GeV. 

If a narrow resonance really exists, there should exist 
(1) a spectral density p(s) with reasonably good behavior 
and (2) some values of So, on which the dependence of me 
is weak. For example, an oscillating p(s) from negative 
to positive values around is regarded as having ’’bad” 
behavior. The physical spectral density should take the 
Breit-Wigner form around m@ if 0 + really exists as a 
very narrow resonance. Hence, p(s) should be either fully 
positive-definite or negative-definite around mg,. 

We will use FESR to analyze the possible existence of 
J = |, S = +1 pentaquarks in this work. This paper is 
organized as follows. In Section El we construct the in¬ 
terpolating currents with different isospin and present the 
formalism. The spectral densities and numerical analy¬ 
sis are given in Section eh The last section is a short 
discussion. 


II. FORMALISM OF FESR 

The method of QCD sum rules incorpo¬ 

rates two basic properties of QCD in the low energy 
domain: confinement and approximate chiral symmetry 
and its spontaneous breaking. One considers a corre¬ 
lation function of some specific interpolating currents 
with the proper quantum numbers and calculates the 
correlator perturbatively starting from high energy re¬ 
gion. Then the resonance region is approached where 


non-perturbative corrections in terms of various conden¬ 
sates gradually become important. Using the operator 
product expansion, the spectral density of the correla¬ 
tor at the quark gluon level can be obtained in QCD. 
On the other hand, the spectral density can be expressed 
in term of physical observables like masses, decay con¬ 
stants, coupling constants etc at the hadron level. With 
the assumption of quark hadron duality these two spec¬ 
tral densities can be related to each other. In this way 
one can extract hadron masses etc. 

We use the following interpolating current for the I = 
0, S = +1, J = | pentaquark state 


V®,k-n( x ) ~ ~^ £ 


Ua(x)Cy 5 d b (x)^ jzi e (x)s e (x)x 
7 M d c (x) - d e (x)s e (x)'y^u c (x ) j (1) 


where three quarks and the remaining qq pair are both 
in a color adjoint representation G2- Similarly, we can 
introduce the I = 1, J = | current 


Vij,,k-n ( x ) — 


Ua(x)Cy 5 d b (x)^ {u e (x)s e (x) 
7 M d c (x) + d e (x)s e (x)7 M u c (x)| . (2) 


For I = 2, I z = 2, J = | state, we use 


vI,k-n ( x ) = \uT(x)Cy„u b {x)\y l ''y 5 u e {x)s e (x)y IJ ,u c (x) . 

(3) 

The overlapping amplitude fj of the interpolating cur¬ 
rent is defined as 


(°M0)|-,p,j) = fjV^ip) 


(4) 


where j is the isospin, v M is the Rarita-Schwinger spinor 
for the J = | pentaquark, which satisfies (p — = 

0, v^ = -2M.v, and 7^ = p^ = 0. 

We consider the following correlation function 


i J d+xe-^ < 0|T{^(x)?7T(0)}|0 >= 

9 y.v (n A (p 2 )p + hl B {p 2 )) H- (5) 

where the ellipse denotes other Lorentz structures which 
receive contributions from both J = 4 and J = | reso¬ 
nances. The tensor structures g flL ,. g IJV p are particular. 
They receive contribution only from the J = § pen¬ 
taquarks. 

We can write a dispersion relation for the scalar func¬ 
tions n A ’ B {p 2 ). 

H a,b{p 2 ) = f ds PA,B ^ S \ (6) 

J (s-p 2 - le) 

where pa,b(s) is the spectral density. 
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The width of the © + pentaquark is less than several 
MeV. For such a narrow resonance, its spectral den¬ 
sity can be approximated by a delta function very well. 
Hence, at the hadron level we have 

Pa{s) = /aH s ~ Mj) + higher states (7) 
Pb(s) = f%Mj5{s — Mj ) + higher states 

where Mj is the pentaquark mass. In principle, there also 
exists non-resonant kaon nucleon continuum contribution 
to Pa,b{s) since 0 + lies above threshold. However, this 
kind of non-resonant K N continuum is either of D-wave 
for J p = | or of P-wave for J p = ^ . Their contribu¬ 
tion is strongly suppressed compared to the resonant 0 + 
pole contribution. Here we want to emphasize that 

fl = +/1 (8) 

for J p = | pentaquarks while 

fl = -fl (9) 


for J p = | pentaquarks. Hence the relative sign be¬ 
tween fg and \Ja\ 2 indicates the parity of the correspond¬ 
ing pentaquark. 

On the other hand, the spectral density pa,b(s ) can be 
calculated at the quark gluon level. For example, the cor¬ 
relation function for the interpolating current 0 reads 


< 0\T{r,l(x)fjt(0)}\0 >= 

^abc^a'b' c' 


{2iS™' C 75 S^ b ' 7 p S^(-x) 7 p S PbP C l5 S a u 

?Tbb' / 


% T nTbc' / 


-is™ c l5 si bb C 75 ^ c luS e s e {-x) lti ST 
-iS e u c ' lv S! e (-xh pSf C 7 5 sl aa ' C 7 5 s b / 


-MTr 


s? aa 'c 75 s b / c 75 s? 7v sf e 

St' e (-x) 7li Sf^ 


->bb', 


-iS e /C 75 Sj bb ' C l5 sf Tr 


-MTr 
—iTr 


s Paa 'c 75 s bb 'c 75 


Tr 




S Paa 'C 75 S b / 7v Sf e (-x) 7 ,SfC 75 Sf } 


-ibc' 


cb', 


where S(x) = — *(0|T{g(x)g(0)}|0) is the full quark prop¬ 
agator in the coordinate space. Throughout our calcu¬ 
lation, we assume the up and down quarks are massless. 
The first few terms of the quark propagator is 


iS ab (x) 


iS a b „ 
2^4 X + 


\n i 

-^gsG^-zia^x + xa^) 


32t r 2 

5 ab S ab x 2 _ 

-u {qq) + -m {9sqaGq) + ' 


First the correlator is calculated in the coordinate 
space. Then n a,b(P) can be derived after making 
Fourier transformation to VLa,b{x). From the imaginary 
part of n A , B (jP) one can extract the spectral density 
Pa,b{s) at the quark hadron level. 


With the spectral density, the nth moment of FESR is 
defined as 


W A ,B(n, s 0 ) 


dss n p A ,B{s ) 


( 10 ) 


where n > 0. With the quark hadron duality assumption 
we get the finite energy sum rule 


WA,B(n, Sq) \Hadron = WA,B(n, So)|qCD • (H) 

The mass and f 2 can be obtained as 


M 2 = W A ,B{n + l,s 0 ) 
4 WA,B(n,s 0 ) 


( 12 ) 


f*Mj = W B (n = 0,s 0 )\ Q cD ■ (13) 

In principle, one can extract the threshold from the re¬ 
quirement 


or 


dM[ 

ds 0 


= 0, 


(14) 


[ (s 0 - s)s n p{s) = 0 . 
Jml 


(15) 


If p(s) > 0 or p(s) < 0 in the whole region [to 2 , oo), there 
does not exist a stable threshold for this finite energy 
sum rule . 

For the gluon condensates we keep only D=4 term and 
neglect D=6, 8 pieces in our calculation. The contribu¬ 
tion of D=6, 8 gluon condensates was found to be much 
smaller than D=4 term in previous QSR analysis. For 
D=7-9 power corrections, we keep only those numeri¬ 
cally large terms, which are related to the quark conden¬ 
sate (qq) or the quark gluon mixed condensate ( g s qaGq). 
Condensates such as (qq)(g 2 GG), (g s qaGq)(g 2 GG ) are 
neglected. 

We use the following values of conden¬ 
sates in the numerical analysis: (qq) = 

— (0.24GeV) 3 , (ss) = -(0.8 ± 0.1)(0.24G eV) 3 , (g 2 GG) = 
(0.48 ± 0.14)GeV 4 , (g s qcrGq) = - nn\ x (qq), 
ml = (0.8 ± 0.2)GeV 2 . We use m s (lGeV) = 0.15GeV 
for the strange quark mass in the MS scheme. 

















III. NUMERICAL ANALYSIS 
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A. 7 = 0 FESR from g Structure 


After tedious calculation, the spectral density p B (s) 
with condensates up to dimension 9 reads 

1 


= 2i».175-X 


+ 25 <s,)) 


7 


2 20 • 45 • 7T 8 
17 


s 3 m s (g 2 GG) 
s 3 (qa ■ Gq ) 


2 14 • 45 • 7r 6 ' 

- 2 12 . 3,^4 s 2 m a(^(gg) - 5 (ss))(qq) 

2 8 . ^7 ^2 s(56(qq) + 75 (ss))(qq) 2 (16) 


where we have used the factorization approximation for 
the high-dimension quark condensates. 

Both the perturbative piece and the D=3 power cor¬ 
rection from the quark condensate in p° B (s) are positive. 
There are two types of D=5 power corrections. One arises 
from the gluon condensate m s (g 2 GG). Its contribution 
is positive. The other one is from the quark gluon mixed 
condensate ( g s qaGq ), which overwhelms the D=3 quark 
condensate and D=4 gluon condensate in magnitude and 
carries a minus sign. As can be seen from Fig. 1, the 
quark gluon mixed condensate renders the spectral den¬ 
sity negative for a big range of s. The D=7 condensate 
yields a positive contribution and cancels the big neg¬ 
ative contribution from the quark gluon mixed conden¬ 
sate, leading to a nearly vanishing p° B {s) for s < 4 GeV 2 
to this order. The contribution from the D=9 condensate 
is also large. In fact, p B (s) > 0 is positive throughout 
the whole range [m 2 ,oo) with the inclusion of the D=9 
power correction. 

It is a common feature in the pentaquark sum rules 
that the quark gluon mixed condensate plays a very strik¬ 
ing, sometimes dominant role when the sum rule is trun¬ 
cated at low orders. In contrast, the ratio between the 
power corrections from ( g s qaGq ) and (qq) is less than 
—20% in the nucleon mass sum rule from p structure j34j. 
which ensures the convergence of the operator product 
expansion (OPE). 

In order to analyze the OPE convergence in the present 
case, we list the right hand side (RHS) of Eq. (Hill for 
the case of n = 0 below: 


W£=°(0, s 0 ) ~ 3.2 x 10 _3 s 0 + 1 - — + H + ^ (17) 

So S 0 s 0 


where the individual terms are normalized according to 
the quark condensate. 

In the FESR framework the expansion parameter is 
1/so while it’s 1/Mg in the conventional Laplace sum 
rule analysis. Eq. GO) indicates that the present FESR 


rho 



FIG. 1: The variation of the spectral density p° B from the 
tensor structure „ (in unit of GeV 11 x 10~ 8 ) with s (in unit 
of GeV 2 ) for the 1=0 current Q. The long-dashed, short- 
dashed and solid curves correspond to p% with D = 5, 7, 9 
condensates respectively. 


is very sensitive to the high dimension condensates with 
D > 5. Numerically, these nonperturbative power correc¬ 
tions converge for s 0 > lOGeV 2 only. Such a big thresh¬ 
old is irrelevant for the 0 + pentaquark around 1.53 GeV. 
We have to conclude that this FESR is not suitable for 
the extraction of pentaquark mass. 

If we ignore the convergence problem, the variation 
of the pentaquark mass with the continuum threshold is 
shown in Fig 2. Naively, there exists a ’’stable” threshold 
of so = 5.6 GeV 2 corresponding to me = (1.8±0.1) GeV 
if we truncate the sum rule at D=5 order. 

M 



FIG. 2: The variation of the mass M with the threshold so 
(in unit of GeV 2 ) corresponding to IIb of the 1=0 current 0. 
Same conventions as in Fig. 1. 


However, if we add the D=7 correction, the shape of 
the curve changes dramatically. With the inclusion of 
D=9 condensates, there is no stable continuum thresh¬ 
old at all. The underlying reason of the extreme sensi¬ 
tivity of this FESR to the high dimension condensates 
is the lack of the convergence of the operator product 
expansion. Hence the extracted pentaquark mass and 
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continuum threshold are very unreliable. 


C. 1=1 and 1=2 Cases 


B. 1 = 0 FESR from g^p Structure 

Now let’s move to the chirally even tensor structure 
g^vP- Pa( s ) with condensates up to dimension 10 reads 

- jdfr^ GG > 

+2win?» 3 ”-(S <5s> - 

:s 2 (^)(qq) 2 + 7(qq)(ss)) 


2 11 • 9 • 7r 4 v 20 


193 


TS 2 m s (qa ■ Gq) 


2 18 • 15 • 7T 6 

+ 2 14 ,27 7 r 4 s ( 1021 (gg) + 286 (ss))(qct • Gq) 
+ ^tt4 - ~m s (-{qq) + ^r(ss))(qq) 2 (18) 


2 6 • 9 • 7T 2 


24 


All the above analysis can be extended to / = 1,2 
case. Roughly speaking, the same conclusions hold. The 
spectral densities for the 1=1 case are 


Pb(s) = 


-1 


and 


2 19 • 175 • 7T 8 


s 5 m s 


w^r^ s ‘ m ~i (5a>) 


1 


2 18 • 15 • 7T 8 
1 

" 2 12 -45-7T 6 


1 


2 8 • 27 • 7T 2 


s 3 m s (g 2 GG) 
s 3 (qa ■ Gq) 

Q7 

s 2 m s {—(qq) - (ss))(qq) 


s(2A(qq) + 29 (ss)){qq) 2 (20) 


We divide the zeroth moment from p\(s) by the D=6 
corrections, which arise mainly from the four quark con¬ 
densate. 

W{=°( 0, s 0 ) ~ 3.8 x 10 _4 Sn + 2 x 10- 2 s 0 +1 - — + ^ . 

s 0 sg 

(19) 

Here the D=8 term, which is the product of (qq) and 
(. g s q<rGq ), plays a very important role. The convergence 
of the OPE requires So > 13 GeV 2 , which renders the 
above FESR useless in the extraction of m©. If we ig¬ 
nore the convergence criteria, we may arrive at rather 
misleading results of m@ and sq as shown by the varia¬ 
tion of with so in Fig. 3 when power corrections with 
D=6, 8, 10 are included. 


Pa(s) = 


11 


23 


2 17 • 5! • 5! • 7T 8 2 22 • 405 • t r 8 

1 3 z 11 / x ! 6 , , 

+ 2^ S 27 SS ~ 15 ^ 


s 3 (g 2 s GG) 


1 


97 

r £S 2 {—(qq) 2 + 3{qq)(ss)) 


89 


TS 2 m s (qa ■ Gq) 


2 18 • 15 ■ 7T 6 
+ 2 13 • 27 ■ 7 r 4 + K(ss))(q* ■ Gq) 

m s (-(qq)+ ^r(ss))(qq) 2 . (21) 


2 6 • 9 • 7T 2 


24 


We also have normalized zeroth moments 


M 



FIG. 3: The variation of M with so corresponding to IIa 
of the 1=0 current 0 . The long-dashed, short-dashed and 
solid curves correspond to p° A with D = 6, 8,10 condensates 
respectively. 


Wt\0,s 0 ) ~ —4.4x 10 -2 so + l ——+ -j + — 5 — • ( 22 ) 

so sg sg 


Wj =1 (0,s o ) ~ 7.7x 10 _4 sg + 7.3x 10 _3 s 0 + l — —+ 


5.1 0.3 

s o 

(23) 


So sg 


For 1=2 case, the spectral densities read 
-1 


Pb( s ) ~ 


4 /- 


2 11 • 27 • 7T 6 ' 
1 

2 15 • 9 • 7T 8 
1 

2 10 • 9 • 7T 6 


(qq) 

s 3 m s (g 2 GG) 
s 3 (qa ■ Gq) 


-72~^4 s 2 m s ((qq) - ss)){qq) 


1 


s{2A{qq) + 23 {ss)){qq) 2 , (24) 


216 ■ 7T 2 
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7 r 59 

= 2 12 -5!-5!-ir sS “ 2» ■ 135 ■ jr 8 ' <9 ' G< ^ 

+ 2l5 i 3?» 3m -<M < “ > - t <99>) 


1 


2 6 • 9tt 4 
499 


2 13 • 45 • 7r 6 


S 2 (2(qq) 2 +3(qq){ss)) 
■s 2 m s (qa ■ Gq ) 


1 751 

- 2 8 , 9 , 7 r 4 ' s (^ 2 ~(gg) + 37(ss))(gcr • Gg) 


1 


■iF^ m>Hs) + 12 


11 J ss))(qqf 


Similarly we have 


Wg 2 (0, So) ~ 1 - \ - 2 

So 


22 505 

J + — ■ ( 25 ) 


M 



FIG. 5: The variation of M with so corresponding to IIa of 
the 1=1 current 0. 


W| =2 (0,s o ) ~ 4.2 x 10 _4 So + 1.8x 10 _2 s o + l — —+ . 

s 0 sg 

(26) 

For completeness we also present the variation of too 
with So in Figs. 4-7 although they are unreliable due to 
the lack of convergence of OPE. 

M 



FIG. 4: The variation of M with so corresponding to IIs of 
the 1=1 current 0 . 


IV. DISCUSSION 

In this paper we have constructed the finite energy sum 
rules for the spin 3/2 pentaquarks using the K-N color- 
octet type interpolating currents. Both g and g fu ,p 
tensor structures are unique for spin 3/2 pentaquarks. 
Because of the high dimension of the pentaquark interpo¬ 
lating currents, power corrections from condensates with 
D > 5 unfortunately turn out to be numerically large for 
our three interpolating currents. Especially the quark 
gluon mixed condensate plays a dominant role. 

Numerical analysis indicates the stable region of the 
continuum threshold is around several GeV 2 . But OPE 


M 



FIG. 6: The variation of M with so corresponding to Ub of 
the 1=2 current 0 . 


of both sum rules converges only when the continuum 
threshold is very big, so > 10 GeV 2 . Such a large value 
of So is irrelevant to the experimentally observed 0 + pen¬ 
taquark. In fact, lack of OPE convergence renders these 
FESRs very sensitive to the high dimension power cor¬ 
rections. 

One may wonder whether the above conclusion is the 
artifact of FESR approach only. We have also performed 
the numerical analysis using the Laplace sum rule. Re¬ 
quiring the pole contribution is greater than 40% of the 
whole sum rule, we arrive at the lower limit of the Borel 
parameter M m ^ n . The convergence of operator prod¬ 
uct expansion requires the high dimension operators be 
suppressed. Numerically we may require the ratio be¬ 
tween the dimension D condensate and perturbative term 
is smaller than \/2 D . In this way we get the upper limit 
of the Borel paramter M m ax- For all the above interpo¬ 
lating currents, we find M m j n > M m ax for both tensor 
structures. In other words, there does not exist a work¬ 
ing Borel window in the Laplace sum rule analysis. This 
point has been noted for the spin 1/2 pentaquark case in 
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M 



FIG. 7: The variation of M with so corresponding to IIa of 
the 1=2 current 0. 

Recall both <f> and to mesons are narrow resonances 
above threshold. 4> decays into KI\ and u> decays into 
37 t. Their FESRs converge. The extracted vector meson 
masses agree with the experimental data. The resonance 
pole contribution dominates the KK or 37r continuum 
(background). Similarly, if spin 3/2 pentaquarks really 
exists as an extremely narrow resonance as indicated by 
those positive experiments, its pole contribution should 
dominate the KN background. Hence one would expect 
(1) a converging FESR at least for one of the two tensor 
structures; (2) a strong signal in the working window with 
the continuum threshold slightly above itiq. But none of 
these spin 3/2 pentaquark FESRs satisfies these condi¬ 


tions. This fact strongly indicates the possible nonexis¬ 
tence of spin 3/2 pentaquarks around 1.53 GeVj which is 
compatible with the most recent CLAS data j23|. With 
a ten times larger database, the nK + spectrum is very 
smooth around 1.53 GeV. In fact, CLAS found no signal 
of exotic baryon resonances with R = +l,S' = +lupto 
2.2 GeV. 

Although our present investigation indicates the pos¬ 
sible nonexistence of a narrow spin 3/2 pentaquark using 
the kaon-nucleon color-octet interpolating currents from 
QCD finite energy sum rule analysis, this is not a strict 
proof yet. An exhaustive study of other interpolating 
currents in search of excellent OPE convergence is very 
desirable. Only after OPE convergence is established, 
may one be able to judge the existence of pentaquarks 
rigorously from the behavior of QCD sum rules. One im¬ 
portant scheme is to include the coupled channel effects 
by introducing a mixed interpolating current which is a 
combination of several interpolating currents with differ¬ 
ent color structures. The mixing and coupled channel 
effect may help suppress the high dimension condensates 
and stabilize the sum rule. Work along this direction is 
in progress. 
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